Glassy dynamics and dynamical 
heterogeneity in colloids 



Luca Cipelletti 

LCVN, UMR 5587 Universite Montpellier 2 and CNRS, P. Bataillon 34095 Montpellier, 

France. 

Eric R. Weeks 

Department of Physics, Emory University: Mail stop 1131/002/1AB 400 Bowman Dr. Atlanta 

GA 30322-2430 USA. 



OXFORD 

UNIVERSITY PRESS 



Abstract 



Concentrated colloidal suspensions are a well-tested model system which has a glass 
transition. Colloids are suspensions of small solid particles in a liquid, and exhibit 
glassy behavior when the particle concentration is high; the particles are roughly anal- 
ogous to individual molecules in a traditional glass. Because the particle size can 
be large (100 nm - 1000 nm), these samples can be studied with a variety of op- 
tical techniques including microscopy and dynamic light scattering. Here we review 
the phenomena associated with the colloidal glass transition, and in particular dis- 
cuss observations of spatial and temporally heterogeneous dynamics within colloidal 
samples near the glass transition. Although this Chapter focuses primarily on results 
from hard-sphere-like colloidal particles, we also discuss other colloidal systems with 
attractive or soft repulsive interactions. 



0.1 Colloidal hard spheres as a model system for the glass 
transition 

0.1.1 The hard sphere colloidal glass transition 

When some materials are rapidly cooled, they form an amorphous solid known as a 
glass. This transition to a disordered solid is the glass transition flGotze and Sjogren, 1992| 
|Stillinger, 199 5; Edig cr et ai, 1996||Angell et ai, 2000] ) . As the temperature of a molec- 
ular glass-forming material is decreased the viscosity rises smoothly but rapidly, with 
little apparent change in the microscopic structure ( Ernst et al., 1991] |van Blaaderen and Wiltzius, 1995 1 
Glass formation may result from dense regions of well-packed molecules or a decreasing 
probability of finding mobile regions. As no structural mechanisms for this transition 
have been found, many explanations rely on dynamic mechanisms. Some theoretical 
explanations focus on the idea of dynamical heterogeneities QGotze and Sjogren, 1992 



Sillescu, 1999 Ediger, 2000 Adam and Gibbs, 19651. The underlying concept is that 



for any molecule to move, all molecules within a surrounding region must "cooperate" 
in their movement. As the glass transition is approached the sizes of these regions 
grow, causing the rise in macroscopic viscosity ( |Adam and Gibbs, 1965[ ). The micro- 
scopic length scale characterizing the size of these regions could potentially diverge, 
helping explain the macroscopic viscosity divergence. However, it is also possible that 
these regions could grow but not be directly connected to the viscosity divergence. 
Additionally, it is not completely clear if the viscosity itself diverges or simply be- 
comes too large to measure ( |Hecksher et ai, 2008[ ). While the existence of dynamical 
heterogeneities in glassy systems has been confirmed in a wide variety of systems, 



the details of this conceptual picture remain in debate (Sillescu, 1999 Ediger, 2000 
Glotzer, 2000} |Ngai, 19991 |Richert, 2002} |Cipelletti and Ramos, 2005[ ) 



Colloidal suspensions are composed of microscopic-sized solid particles in a liquid, 
and arc a useful model system for studying the glass transition. In terms of intcrparticlc 
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interaction, the simplest colloids are those in which the particles interact as hard 
spheres, i.e., the interparticle potential arises solely due to excluded volume effects 
(Pusey and van Megen, 19861. Hard spheres are a useful theoretical model for glass- 



forming systems due to their simplicity (Bernal, 1964). Clearly, attractive interactions 
between atoms and molecules are responsible for dense phases of matter. But given 
dense states of matter, repulsive interactions play the dominant role in determining 
the structure. Hard spheres are useful simulation models for crystals, liquids, and 
glasses, although it is still debated whether a purely repulsive interparticle potential 



is sufficient to reproduce the glass transition in general ( |Berthier and Tarjus, 2009 1 . 

The control parameter for hard sphere systems is the concentration, expressed as 
the fraction, ip, of the sample volume occupied by the particles. Most colloidal hard 
sphere systems act like a glass for ip larger than an "operational" glass transition 
volume fraction ip g « 0.58. The transition is the point where particles no longer dif- 
fuse through the sample on experimentally accessible time scales; for p < <p g spheres 
do diffuse at long times, although the asymptotic diffusion coefficient decreases 
sharply as the concentration increases dBartsch et al, 1993[ | van Megen et al, 1998[ 
|Kasper et al, 199"8| . The transition at ip g occurs even though the spheres are not 
completely packed together; in fact, the density must be increased to </?rcp ~ 0.64 for 



"random-close-packed" spheres ( O 'Hern et al, 2003 Bernal, 1964 Torquato et al, 2000 



Donev et al, 2004 O'Hern et al, 2004) before the spheres are motionless. In simula- 
tions, a collection of same-sized (i.e. monodisperse) hard spheres almost always crystal- 
lizes. A binary mixture of spheres or indeed a distribution of sizes is therefore needed 
to frustrate crystallization and enable access to the glass transition, both numerically 
and experimentally ( |Henderson and van Megen, 1998 Zaccarelli et al, 2009a). 



van Megen and Puse 


y, 1991 


van Blaaderen and Wiltzius, 1995||van Megen and Underwood, 1993[ 


Bartsch et al., 1993 


Bartsc 


l, 1995 Mason and Weitz, 1995) is comparable in many 



Macroscopically, a colloidal liquid flows like a viscous fluid whereas the colloidal glass 
does not flow easily, like a paste ( |Segre et al, 19956; Chen g et al., 2002| . For colloidal 



samples with low polydispersity, samples can crystallize for ip < ip g (with crystalliza- 
tion nucleating in the interior of the sample), while for cp > tp g , crystals only nucleate at 
flat sample boundaries such as the walls of the container flPusey and van Megen, 19861 . 
Microscopically, the glass transition point is identified as the point where — > 



(Bartsch et al., 1993 van Megen et al., 1998 Kasper et al., 1998) 



Note that there are some questions about the colloidal glass transition. First, 
prior measurements disagree about the nature of the viscosity divergence in col- 
loidal glasses QSegre et al, 1995~6| |Cheng et al., 20 02). This may be due to difficul- 
ties in reconciling measurements of the volume fraction ip ( |de Schepper et al., 1996{ 
Segre et al., 1996[) . Second, it has been seen that a glassy colloidal suspension can crys- 



tallize under microgravity conditions (Zhu et al., 1997), and potentially also when in 
density-matching solvents (Ke gel, 2000[ " Simeonova and Kegel, 2004[ ), suggesting that 
the apparent glass transition at ip g = 0.58 is an artifact of gravity. The interpretation 
of these observations is unclear. One possibility is that this is merely heterogeneous 
nucleation at the walls, as seen before ( |Pusey and van Megen, 1986[ ). Another possi- 
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bility is that these samples were more monodisperse than most, and a slightly larger 
polydispersity (> 5%) may be necessary to induce a glass transition at <p g = 0.58 
(|Henderson and van Megen, 199 8 ; Mell er and Stavans, 1992[ |Auer and Frenkel, 2001 



|Schope et ai, 2007 ). In practice, most experimental samples always have a polydisper- 
sity of at least 5%. Recent simulations suggest that the relationship between polydis- 
persity, crystallization, and glassy dynamics is more complex than perhaps was pre- 
viously appreciated ( |Zaccarelli et ai, 2009a[ ) . Overall we note that the interpretations 
of the various observations described in this paragraph are often still controversial. 
Colloidal glasses are quite similar to molecular glasses: 



Both are microscopically disordered ( |van Blaaderen an d Wiltz ius, 1995 1 



Both are macroscopically extremely viscous near their transition point (Scgre et ai, 1995 & 



Cheng et ai, 2002 1 



Colloidal glasses have a nonzero elastic modulus at zero frequency, which is absent 
in the liquid phase ( |Mason and Weitz, 19951 ). 

Colloidal glasses are out of equilibrium and show aging behavior - their proper- 
ties depend on the time elapsed since preparation (Courtland and Weeks, 2003{ 



Cianci et ai, 2006 Cianci and Weeks, 2007 Simeonova and Kegel, 2004 1 , similar 



to polymer glasses ( Hodge, 1995 |McKenna, 2003[ ) and other molecular glasses 



dAngell et ai, 2000| |Castillo and Parsaeian, 2007 



Colloidal glasses exhibit dynamical heterogeneity ( Kegel and van Blaaderen, 2000 



Weeks et al., 2000 1 , similar to that seen in simulations (Donati et ai, 1998: Glotzer, 2000 



Doliwa and Heuer, 1998) and experiments on molecular glasses (Ediger, 2000 



Sillescu, 1999] [Vidal Russell and Israeloff, 2000 1 



• The colloidal glass transition is sensitive to finite size effects ( |Nugent et al., 2007| ), 

similar to molecular glasses ( |Alcoutlabi and McKcnna, 2005 ) and polymers ( Roth and Dutcher, 2005[ ) . 

One difference between colloids and molecules is that colloidal particles move via 
Brownian motion whereas the latter move ballistically at very short time scales; sev- 
eral simulations indicate that this difference is unimportant for the long-time dy- 
namics which is of the most interest ( |Gleim et ai, 1 998; Sz amel and Flenner, 2004"[ 
|Hofiing et al., 2008[ ). Likewise, hydrodynamic interactions between particles influence 
their motion on short time scales, but do not modify the pairwise interaction potential 
(which remains hard-sphere- like), suggesting that they should not be relevant for the 
long-time dynamics {Brambill a et ai, 2009[ ). 

One advantage colloids have over traditional molecular glassformers is that their 
time scales are significantly slower, with relaxation taking 0(1-1000 s), allowing easy 
study of the relaxation processes. A second advantage of colloids is that their large 
size [0(1 /xm)] allows for measurement using optical microscopy or dynamic light 
scattering, as will be discussed in Sec. 10.21 

Experiments have one chief advantage over simulations, in that they more easily 
avoid finite size effects. Near the glass transition, dynamical length scales can be 
large (^4—5 particle diameters ( [Weeks et al., 2 007 Do liwa and Heuer, 2000[ )) and 
finite size effects on structure and dynamics may extend to even larger length scales 
(~ 20 particle diameters ( Nugent et ai, 2007[ )). Microscope sample chambers typically 
contain ~ 10 9 particles, and light scattering cuvettes contain even more. 

The most popular "hard sphere" colloid is colloidal PMMA (poly- methyl- mcthacrylatc) , 
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sterically stabilized to minimize inter-particle attraction flPusey and van Megen, 1986| 



| Weeks et al., 2000] |Dinsmore et al., 2001| |Antl et al., 1986 1 . The particles can be placed 



in density-matching solvents to inhibit sedimentation (see Sec. 10.1. 2"]) . and/or solvents 
which match their index of refraction to enable microscopy ( [Dinsmore et al., 2001 ) 



and light scattering (see Sees. 10.2^21 and [0~.2.3[) . In some of these solvents, the particles 
pick up a slight charge and thus have a slightly soft repulsive interaction in addition to 
the hard-sphere core. Despite this charge the spheres behave similarly to hard spheres 
with some phase transitions shifted to slightly lower ip ( |Gasser et al., 200~lj ). Salt can 
be added to the samples to screen the charges and shift the interaction back to more 
hard-sphere like ( |Yethiraj and van"Fj laadcrcn, 2003; Royall et al., 2003). 



A second popular colloid is colloidal silica, which is relatively easy to fabricate 
( |Stober, "19 68). These particles are suspended in water (or a mixture of water and 
glycerol), avoiding the organic solvents that are required for PMMA colloids. Be- 
cause they are in water, repulsion due to charge is the primary mechanism preventing 
flocculation; adding salt can screen the charges and cause flocculation (which is often 
irreversible). Silica colloids are hard to density match (their density varies but is larger 
than 2 g/cm 3 ), and also it is hard to match their refractive index. This latter constraint 
makes microscopy difficult except at lower volume fractions ( |Mohraz et al., 2008[ ). 

Several glass transition theories have been applied to the colloidal glass transi- 
tion. The colloidal glass transition appears to be well-described by mode-coupling the- 



ory up to ip « 0.58 ( van Megen and Pusey, 1991 van Megen and Underwood, 1993 
van Megen and Underwood, 1994||Gotze and Sjogren, 1991||Schweizer and Saltzman, 2003 



Schweizer and Saltzman, 2004; Saltzman and Schweizcr, 20066; Saltzman and Schwci zer, 2006 a| 



although other glass transition theories successfully capture many features of the 
colloidal glass transition as well flNgai and Rendell, 1998| |Liu and Oppenheim, 1996[ 
Tokuyama and Oppenheim, 1995[ Tokuyama, 2007[ ). (See the discussion in Sec. 10.3.11 



which discusses the strengths and weaknesses of mode-coupling theory as applied to 
colloids.) 



0.1.2 Experimental challenges in studying colloidal hard spheres 

As argued in the previous section, colloidal hard spheres are a good model system 
for investigating the glass transition. However, several experimental challenges have 
to be faced. Probably, the most serious problem is that of a precise determination 
of the volume fraction. In optical microscopy, tp can be obtained by counting the 
number of particles in a given volume and using the particle size as obtained, e.g., 
from electron microscopy. It should be noted that an error of just 1% in the radius a 
of the particles results in a 3% error in ip. Additionally, one has to take into account 
the thickness of the stabilizing layer (e.g. the grafted polymer for PMMA particles 
or the counterion cloud for silica particles in a polar solvent), which is often difficult 
to measure precisely. Since typical values of the thickness are on the order of 10 
nm ( |Pusey, 199ip , this contribution can be relevant for the small particles used in 
light scattering (a ~ 100 — 300 nm), while it is less important for the micron-sized 
particles used in microscopy. Other methods to determine tp include precise density or 
refractive index measurements ( |Phan et al., 1996| . For PMMA, these methods require 
special care, since the particles can be swollen by organic solvents such as tetralin or 
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brominated solvents, which change their density and refractive index compared to 
those of the bulk material. 

Because of these difficulties, the absolute volume fraction is often determined in- 
directly by comparing the phase behavior or the dynamical behavior of the sam- 
ple to theoretical and numerical predictions. Samples sufficiently monodispcrsc [a = 

Va?-a 2 /a 2 < 5-8%) crystallize for 0.494 < p < 0.545 QPusey, 199ft . Th e absolute 
volume fraction can then be calibrated by matching the experimentally determined 
freezing volume fraction with ip f 



0.494 as determined by simulations (Hoover and Ree, 1968 1 



However, some uncertainty is still left, because the exact value of the freezing fraction 

depends on a ( |Bolhuis and Kofke, 19 96 Fasol o and Sollich, 2004|[van Megen and Underwood, 1994[ 



|Zaccarelli et ai, 2009a| Pusey et al., 2009 ). Moreover, this method cannot be applied 



to more polydisperse suspensions that do not crystallize over several months or years. 
Alternatively, p may be calibrated against predictions for the volume fraction depen- 



dence of the low shear viscosity ( Pusey, 1991| Poon et al., 19961 or the short time self 
diffusion coefficient ( |Beenakker and Mazur, 1983| Tokuyama and Oppenheim, 1994[ ) 
in the dilute regime. For samples where both the calibration against ipf and that 
using the short time self diffusion coefficient are possible, the two methods appear 
to be consistent (van Megen and Underwood, 1989; Segre et al., 1995a|. In summary, 



while relative values of p can be measured very precisely (down to 10 using an an- 
alytical balance), absolute values are typically affected by an uncertainty of a few %. 
This should always be kept in mind when comparing sets of data obtained in different 
experiments. 

For colloids, the equivalent of T = in a molecular system is random close packing, 
the volume fraction </Jrcp ~ 0.64 where osmotic pressure diverges and all motion ceases 
because no free volume is left. Clearly, knowledge of the precise location of </?rcp is 
very important to discriminate between theories that predict a glass transition below 



random close packing (for example the mode coupling theory (Gotze, 1999) or thermo- 
dynamic glass transition theories ( |Cardenas et aL,T9 99; Parisi and Zamp om, 2005[ )) 
and scenarios like jamming, where no arrest is predicted below pnep- However, the 



location of random close packing is still highly debated ( Berthier and Witten, 2009 



|Xu et al, 2009; K amien and Liu, 2007| |Chaudhuri et al., 2010 1 and its very existence 
is challenged, based on the argument that one can always trade order for packing 
efficiency flDonev et al., 2007] ), up to p = n/y/l8 ps 0.7405, the packing fraction 
of a (monodisperse) hard sphere crystal. Experimentally, it is difficult to measure 
Vrcp because of the uncertainty on absolute volume fractions discussed above and 
because applying the high pressure (e.g. by centrifugation) needed to approach it 
may result in the compression of the stabilizing layer. Additionally, both experi- 
ments ( |Clusel et al, 2009[ ) and simulations ( |Schaertl and Sillescu, 1994| [Berthier and Witten, 2009| ) 
show that tpRCP depends on a. 

Some of the experimental challenges posed by hard spheres stem from the very 
same features that make them a valuable model system: their relatively large time 
and length scales. The microscopic time in a colloidal system is the Brownian time 
tb , i.e. the time required by a particle to diffuse over its own size in a diluted system, 
defined as 



tb — a 2 /6D = mja 3 /ksT . 



(0.1) 
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Here, rj is the solvent viscosity, T the absolute temperature, ks Boltzman's constant 
and D is the diffusion coefficient for a sphere of radius a, given by the Stokcs-Einstcin- 
Suthcrland formula ( |Einstein, 19 05; Sutherland, 19051: 



D = k B T/6irr)a . (0.2) 

The time scale tb is typically on the order of 10 -3 — Is. The largest relaxation time 
that can be measured in concentrated systems is on the order of 10 5 s; thus, the acces- 
sible dynamical range covers at most 8 decades, as opposed to 15 decades in molecular 
glass formers. Additionally, owing to their relatively large size and due to the mis- 
match between their density and that of the solvent in which they are suspended, 
colloidal particles experience gravitational forces that can modify their phase behav- 
ior ( |Zhu et ai, 1997||Pusey et ai, 2009[ ) and dynamical properties ( |E1 Masri et ai, 2009} 



|Simeonova and Kegel, 20*04 1 . The relevant parameter to gauge the importance of sed- 



imentation is the inverse Peclet number, Pe 1 = , 3 . fcBT 4 , defined as the ratio of the 
gravitational length to the particle radius (g is the acceleration of gravity and Ap 
the density mismatch). For diluted suspensions, gravity becomes relevant as Pe -1 ap- 
proaches (from above) unity; for concentrated suspensions, sedimentation effects may 
set in at even higher values of Pe -1 , since gravitational stress is transmitted over 
increasingly larger distances as the sample becomes more solid-like. For example, sed- 
imentation effects have been reported to alter the dynamics of PMMA particles in 
organic solvents (Ap ~ 0.3g/cm 3 ) for Pe -1 = 44.1 (S imeonova and Kegel, 2004] ) or 



even for Pe = 1350 (El Masri et ai, 20091. Density- matching solvents can miti 



gate these effects, although matching closely both the index of refraction (as required 
for optical observations) and the density of the particles without altering their hard 
sphere behavior has been proved difficult for PMMA particles flRoyall et ai, 2003] ) and 
impossible for other systems such as silica spheres. 

0.2 Experimental methods for measuring both the average 
dynamics and dynamical heterogeneity 

0.2.1 Main features of optical microscopy and dynamic light scattering 

Optical microscopy and dynamic light scattering are the main techniques to probe 
both the average dynamics and its spatiotemporal fluctuations in dense colloidal sus- 
pensions. Each of them comes with specific advantages and limitations. Optical mi- 
croscopy is unsurpassed in providing detailed information on the structure and the 
dynamics at the single particle level. The same quantities introduced in theory and 
simulations to characterize the dynamics can be precisely measured (e.g. the mean 
square displacement or the intermediate scattering function for the average dynamics, 
and the dynamical susceptibility, %4i an d the spatial correlation of the dynamics, g±, 
for its fluctuations, see Chapter 2 and Sec. 10.2.41 below). Additionally, direct visual- 
ization of the sample allows any experimental problem to be readily detected, such as 
particle aggregation, sedimentation, or wall effects. Finally, techniques such as opti- 
cal or magnetic tweezing QGrier, 2003| Amblard et ai, 1996[ ) allow one to manipulate 



single particles and thus to measure the microscopic response of the system to a local 
perturbation ( |Habdas et ai, 2004] ) . 
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Dynamic light scattering probes a very large number of particles simultaneously, 
yielding very good averages. Moreover, particles used in light scattering are usually 
smaller than those for optical microscopy (a = 100 — 500 nm as opposed to a = 
0.5 — 1.5/im), which has a twofold advantage. First, the microscopic time tb (see 
Eq. (|0.1jl ) is significantly reduced, since tb ~ a 2 / D ~ a 3 , thus increasing substantially 
the experimentally accessible dynamical range. Second, gravitational effects are much 
less of concern, since Pe^ 1 ~ a~ 4 . Finally, as we will discuss it in Sec. 10.2.41 recent 
developments allow dynamical heterogeneity to be probed by dynamic light scattering, 
although not at the level of microscopic detail afforded by optical microscopy. These 
methods extend the possibilities of light scattering by adding features characteristic 
of imaging techniques. Quite in a symmetric way, very recent microscopy methods 
such as dynamic differential microscopy ( |Cerbino and Trappe, 2008[ ) have extended 
imaging techniques by adding the capability of measuring the intermediate scattering 
function f(q, r) defined below in Sec. 10.2.31 

0.2.2 Optical and confocal microscopy 

Given the large size of many colloidal systems (particle radius a ~ 100 — 1000 nm), 
optical microscopy is a useful tool for observing these systems. First, these sizes are 
comparable to the wavelength of light, thus rendering them visible. Second, the time 
scales of their motion are often slow enough for video cameras to follow their motion. 
This can be seen by considering how quickly particles diffuse. Colloidal particles un- 
dergo Brownian motion due to thermal fluctuations, as discussed in Sec. 10.1.21 The 
diffusion coefficent D (Eq. I0.2[) is related to the mean square displacement of particles 
as: 

(Ax 2 ) = 2D At, (0.3) 

where At is the time scale over which the displacements are taken. A particle of diam- 
eter 2a = 1 fjxci in water (77 = 1 mPa-s) at room temperature diffuses approximately 
1 /xm in 1 s. This motion is easy to study with a conventional video camera and a 
microscope; most video cameras take data at 30 images per second, thus they allow 
one to follow the Brownian motion of colloidal particles of these sizes. Of course, par- 
ticles that are 10 times smaller move 1000 times faster, by Eq. (10. ip . In practice, often 



one can choose to study larger colloidal particles (Weeks et al., 20001 or use careful 
data analysis techniques to learn information about smaller-sized particles that aren't 
directly imaged ( |Simeonova and Kegel, 2004| . Additionally, when studying dense col- 
loidal samples, the time scales increase simply because of the glassy dynamics. 

Two main methods of microscopy have been used to study dense colloidal samples: 
conventional optical microscopy and confocal microscopy. First, there is the possibility 
of using a conventional light microscope technique such as "brightfield microscopy." 
These techniques typically depend on slight differences between the index of refraction 
of the colloidal particles and the solvent ( jlnoue an d Sprin g7T-997[ ) . A limitation is that 
these differences also scatter light; each particle acts like a tiny lens. This ultimately 
limits how deeply into a sample one can observe. (This is the same phenomenon 
that makes milk appear white, even though composed of transparent components; 
the different components all have different indices of refraction. Snow is white for 
a similar reason, due to the contrast in index of refraction between the ice crystals 
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and air.) A further limitation of conventional optical microscopy is that the images 
are limited to a plane, although this can be fine for quasi-two-dimensional samples 
( |Marcus et al, 19990 . 

A related conventional technique is fluorescence microscopy. Here, the particles 
can be precisely index-matched with the solvent. However, the particles also con- 
tain a fluorescent dye. In fluorescence microscopy, the particles are illuminated with 
short- wavelength light. The dye molecules absorb this light, and radiate slightly longer 
wavelength (lower energy) light, which is imaged by the camera. Special filters and 
mirrors are used to direct the light appropriately from the light source to the sam- 
ple, and from the sample to the camera. While this method avoids the problem of 
light scattering off of different parts of the sample, in dense samples it can still be a 
problem that too much of the sample fluoresces at the same time, thus giving a large 
background illumination. Trying to observe bright particles on a bright background 
thus limits fluorescence microscopy of dense samples. One way to overcome this is to 
only dye a few tracer particles. 

Fluorescence microscopy has one significant limitation: photoblcaching. After dye 
molecules absorb the excitation light, but before they emit light, they can chemically 
react with oxygen present in the sample to form a non-fluorescent molecule. This only 
happens when they are excited, so photobleaching happens in direct proportion to the 
illumination light. Photobleaching manifests itself as the image becoming gradually 
darker. This can be a useful technique for studying local diffusion in samples, a tech- 
nique known as "fluorescent recovery after photobleaching" ( |Axelrod et al., 1976[ ) . In- 
tense light is used to photobleach a region of the sample, and then low-intensity light is 
used to monitor the recovery of fluorescence as non-bleached particles diffuse back into 
the region. With this method, the diffusivity of the particles can be measured, which 
has been used to study the behavior of colloidal glasses ( Simeonova and Kegel, 2004 ). 

An extension of fluorescence microscopy is confocal microscopy, sometimes termed 
laser scanning optical microscopy. Here, a laser is used to excite fluorescence in dye 
added to a sample. Typically, the laser beam is reflected off two scanning mirrors that 
raster the beam in the x and y directions on the sample. Any resulting fluorescent light 
is sent back through the microscope, and becomes descanned by the same mirrors. A 
mirror directs the fluorescent light onto a detector, usually a photomultiplier tube. 

One additional modification is necessary to make a confocal microscope: before 
reaching the detector, the fluorescent light is focused onto a screen with a pinhole. 
All of the light from the focal point of the microscope passes through the pinhole, 
while any out of focus fluorescent light is blocked by this screen. This spatial filtering 
technique blocks out the background fluorescence light, allowing the particles to be 
viewed as bright objects on a dim background. 

This ability to reject out-of-focus fluorescent light directly results in the main 
strength of confocal microscopy, the ability to take three-dimensional pictures of sam- 
ples. By rejecting out-of-focus light, a crisp two-dimensional image can be obtained, 
as shown in Figure 10.11 The sample (or objective lens) can be moved so as to focus 
at a different height z within the sample, and a new 2D image obtained. By collect- 
ing a stack of 2D images at different heights z, a 3D image is built up. The time 
to scan one 2D image can range from 10 ms to several seconds, depending on the 
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Fig. 0.1 Image of 2 micron diameter colloidal particles taken with a confocal micro- 
scope. The scale bar is 10 microns. The sample is in coexistence between a colloidal crystal 
and a liquid, see Ref. ( |Hernandez-Guzman and Weeks, 2009) for details. Taken by Jessica 
Hernandez-Guzman and Eric R. Weeks 

details of the confocal microscope and the desired image size and quality. The time 
to scan a 3D image depends on the 2D scan speed and the desired number of pixels 
in the z-direction; reasonable 3D images can be acquired in 2 - 20 s depending on 
the microscope dNugent et ql, 2007| [Weeks et ql, 20000. Finally, we mention Coher- 
ent anti-Stokes Raman scattering (CARS) microscopy ( |Kaufman and Weitz, 2006] ), a 
technique that allows one to image in 3D colloidal samples with spatial and temporal 
resolutions comparable to those of confocal microscopy. 

Once a series of images has been acquired (2D or 3D), the next step is typi- 
cally tracking the individual particles within the images. Within each image, a com- 
puter can determine the positions of all of the particles. If the particles do not 
move large distances between subsequent images, then they can be easily tracked 
( |Crocker and Grier, 1996[ ). Specifically, they need to move less between images than 
their typical inter-particle spacing. With confocal microscopy, particles can be tracked 



in three dimensions (Bcsseling et ql., 2009 Dinsmore et ql., 2001). This then is the 



same type of data that is analyzed from computer simulations. Simulations, of course, 
have many advantages, including a tunability of particle interaction, the ability to 
study either Brownian or ballistic dynamics, and more precise and instantaneous con- 
trol over parameters such as temperature and pressure. The experiments have an 
advantage that typically the boundaries are far away: while perhaps a few thousand 
particles might be viewed, they are embedded in a much larger sample with millions 
of particles. 

0.2.3 Dynamic light scattering 



Dynamic light scattering (DLS) ( |Berne and Pecora, 1976 1, also termed photon cor 



relation spectroscopy, probes the temporal fluctuations of the refractive index of a 
sample. In colloidal systems, scattering arises from a mismatch between the index of 
refraction of particles and that of the solvent, so that DLS probes particle density flue- 
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tuations. Experimentally, one measures (72 ("J") — 1, the time autocorrelation function of 
the temporal fluctuations of the intensity scattered at a wave vector q = Air/ A sin(0/2) . 
Here, 9 is the scattering angle and A is the wavelength in the solvent of the incoming 
light, usually a laser beam. Under single scattering conditions, the intensity autocorre- 
lation function is directly related to the intermediate scattering function f(q,r) (ISF, 
sometimes also referred to as the dynamic structure factor) : 



f(q, t) = ( iV- 1 £ exp {-zq ■ [ Tj (t + r) - r k (t)}} 
\ 3,k I 



x //3- 1 [ff2(r)-l] = Jp- 



< I(t + r)I(t) > t 

< m >? 



(0.4) 



where (3 < 1 depends on the collection optics, I(t) is the time varying scattered 
intensity, Vj the position of the j'-th particle, and N the number of particles in the 
scattering volume. Note that f(q, r) decays significantly when Ar(r) = r(i+r) — r(i) is 
of the order of depending on the choice of 9, DLS probes motion on length scales 
ranging from tens of nm to tens of /im. Another important point to be noticed is that 
DLS usually probes collective motion, since the sum in Eq. (|0.4[) extends over all pairs 
of particles. However, there are ways to measure the self part of the ISF by making the 
contribution of the j =/= k terms vanish from Eq. (j0.4p . This may be accomplished by 
choosing the scattering vector q in such a way that S(q) = 1, where S(q) is the static 
structure factor ( |Pusey et ai, 1982] |Pusey, 1978[ ) . Alternatively, one can use optically 



polydisperse suspensions, as in Ref. ( van Megen and Underwood, 1989[ ), where silica 



and PMMA particles of nearly the same radius but different refractive index were 
mixed. When matching the average refractive index of the colloids, the measured ISF 
contains only contributions from the self terms, i = k. Note that optically composite 
particles that are polydisperse in size are usually also optically polydisperse. This is 
the case, e.g., of PMMA colloids stabilized by polymer layer with a refractive index 
different from that of the core, when the core is polydisperse in size. 

The nature of the averages indicated by the brackets in Eq. (|0.4I) is an important 
issue. In the definition of the ISF, the average is over an ensemble of statistically equiv- 
alent particle configurations, while operationally g2 is averaged over time. Therefore, 
ergodicity is required for Eq. (|0 .4[) to hold. Additionally, in order to reduce noise to an 
acceptable level, gi has to be averaged over at least 10 3 — 10 4 r a> with r a the relaxation 
time of the ISF. These requirements are often impossible to meet for supercooled or 
glassy colloidal systems, where r a can be as large as hundreds of thousands of sec- 
onds. To overcome these difficulties, various schemes have been proposed, among which 
the most popular is probably the "multispeckle" method QWong and Wiltzius, 1993| 
|Bartsch et al., 1997[ ). In a multispeckle experiment, the phototube or avalanche pho- 
todiode used in regular DLS is replaced by a multielement detector, typically a CCD 
or CMOS camera sensor. The collection optics is chosen in such a way that each pixel 
of the detector corresponds to a different speckle ( Goodman, 2007[ ), i.e. to a slightly 



different scattering vector. Because distinct speckles carry statistically independent 
information, the time average can be replaced in part by an average over pixels: 
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X < l p (t + t) > p < I p (t) > p 



Here I p indicates the intensity measured by the p — th pixel, < • • • >t and < ■ ■ ■ > p 
denote averages over time and pixels, respectively. The set of pixels is chosen in such 
a way that they correspond to nearly the same magnitude of the scattering vector q. 
The number of pixels is typically of order 10 4 — 10 5 , so that time averaging only needs 
to extend over a few r Q . This approach allows very slow and non-stationary dynamics 
to be probed effectively. 

Although many DLS experiments have been carried on a variety of glassy colloidal 
systems, the single-scattering conditions required by this technique are probably more 
the exception than the rule. For mildly turbid suspensions, smart detection schemes, 
most of which were pioneered by K. Schatzel ( |Schatzel, 1991[ ), allow the rejection of 
multiply scattered photons, thereby efficiently suppressing artifacts due to multiple 
scattering. Popular implementations of this concept include the so-called two-color 
and 3-D apparatuses (see Ref. QPusey, 19 99) for a review). For very turbid samples, 
where photons are scattered a large number of times before leaving the sample and 
the contribution of single scattering is negligible, an alternative formalism has been 
developed, termed Diffusing Wave Spectroscopy (DWS) ( |Weitz and Pine, 1993 ). In a 



DWS experiment the intensity correlation function gi — 1 is related to the mean squared 
displacement, < Ar 2 (r) >, rather than to the ISF, as in DLS. Another important 
difference is the probed length scale, which in DWS typically covers the range 0.1 - 
100 nm, much smaller than in DLS. Finally, DWS experiments typically probe the self 
motion of the particles, rather than their collective relaxation. 

An alternative way to tackle multiple scattering is provided by X-photon correlation 
spectroscopy (XPCS). Modern synchrotron sources deliver X-ray radiation that is 
coherent enough to perform the same kind of experiments as with a laser beam in DLS. 
Because colloidal systems scatter X-ray much less efficiently than visible light, in most 
cases XPCS measurements can be safely performed in the single scattering regime 
without adjusting the refractive index of the solvent. While long-term beam stability 
is often still an issue, several XPCS studies on the slow dynamics of colloidal systems 



have been published in the last years (Bandyopadhyay et al, 2004 


Chung et al, 2006 


Robert et al, 2006||Trappe et al, 2007| [Wandersman et al, 2008 


Herzig et al, 2009 


Duri et al, 2009 a|). 



0.2.4 Time and space resolved dynamic light scattering 

In a traditional DLS experiment, the detector is placed in the far field, so that it 
collects light scattered by a macroscopic region, typically of volume 1 mm 3 or more. 
Additionally, the intensity correlation function 52 ~ 1 has to be extensively averaged 
over time. Because of these averages over both time and space, no information can 
be a priori extracted on dynamical heterogeneity, e.g. on the spatial and temporal 
fluctuations of the dynamics. In recent years, however, novel light scattering methods 
have been proposed to overcome these limitations, providing either spatially averaged 
but temporally resolved data (time resolved correlation, TRC QCipelletti et al, 2003[ )), 
or both spatially and temporally resolved measurements (photon correlation imaging, 
PCI jDuri et al, 2009&1 )). 



xvi Abstract 



In a TRC experiment, one uses a CCD or CMOS detector to calculate a two-time 
correlation function Ci(t, r) defined by 

ej(f>r)= <I P (t + r)I p (t)>» x 

V ' <Ip(t + T) > p <I p (t) > p V ^ 

Note that the usual intensity correlation function <?2(t) — 1 defined in Eq. (|0.5[) is the 
temporal average of ci(t, r). Because the detector is typically placed in the far field, cj 
is a temporally resolved but spatially averaged correlation function. Figure 10.21 shows 
an example of TRC data and their relationship to gi — 1 for a diluted Brownian sus- 
pension (Duri et al., 20051 and for a colloidal gel QDuri and Cipelletti, 2006 ). When 



plotted as a function of time t for a fixed time delay r, the data for the Brownian sus- 
pension are essentially constant. Indeed, for this system the dynamics are homogeneous 
and timc-translational invariant, so that the evolution of the system and hence the de- 
gree of correlation over a fixed time lag r does not depend on t. The small fluctuations 
around the mean value are due to the statistical noise of the measurement, associated 
with the finite number of pixels over which c/ is averaged ( |Duri et al., 2005] ). For the 
colloidal gel, by contrast, c/ exhibits significant temporal fluctuations, indicative of 
heterogeneous dynamics. Sudden drops of cj measured for short lags, as in the top 
trace of Fig |0.2b . arc indicative of a sudden rearrangement event that has led to a loss 
of correlation between the intensity patterns recorded at times t and t+T. As the rear- 
rangement event ceases, the degree of correlation recovers its typical level. At longer 
time lags (middle trace in Fig |0.2b ), cj has a highly fluctuating behavior, because sev- 
eral events may occur during the probed lag. By contrast, almost no fluctuations are 
observed at very long lags (bottom trace in Fig |0.2b ), since a large number of events 
has occurred for all pairs of images, leading to a full decorrelation. 

Various way of characterizing the fluctuations of cj have been proposed, including 
analyzing the probability distribution function and the moments of its fluctuations, or 
their temporal autocorrelation ( |Duri et al., 2005] ). Here, we focus on the variance 

X (t) =var( Cj ) = ([c/(t,r)- < c 7 (t,r) > t )] 2 ) 4 , (0.7) 

which is the analogous in light scattering of the dynamical susceptibility Xi discussed 
in detail in Chapter 2. Intuitively, one understands that large fluctuations of c/ must be 
associated to "rare" , large events: if the rearrangements were very localized, many such 
events would be necessary to significantly decorrelate the light scattered by a macro- 
scopic sample volume. Assuming independent events, the resulting spatial average 
would yield a smooth c/ trace. More precisely, x can be shown to be proportional to the 
volume integral of the spatial correlation of the dynamics, as discussed for Xi in Chap- 
ter 2. An example of the scaling of x with the size of the events for a coarsening foam, 
where events can be unambiguously identified, is discussed in ( |Mayer et at, 2004] ) (see 
also Chapter 5). A few differences exist between x i n light scattering experiments and 
X4 in simulations or real space measurements. Contrary to Xi> X is n °t normalized 
with the respect to the number N of particles (compare Eq. (|0.7[) to eqn (11) in Chap- 
ter 3), since usually N is not known precisely in light scattering. Accordingly, typical 
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Fig. 0.2 a) Degree of correlation ci(t,r) for a diluted suspension of Brownian particles: 
the dynamics are stationary and homogeneous, as seen by the very small fluctuations of cr, 
due uniquely to measurement noise. From top to bottom, r = 0, 1, 2.5, 5, 10 and 700 sec. b) 
Intensity correlation function — 1 obtained by averaging over time the data in a). The 
solid circles correspond to the time delays for which ci is shown in a), c), d): degree of 
correlation and intensity correlation function for a colloidal gel ( |Duri a nd Cipcllctti, 2006) 
(see Sec. 10. 5. 2]) . In c), t = when the gel is formed. Note the large fluctuations of ci, due to 
the heterogeneous nature of the dynamics. Individual events are discernable in the top trace 
(r = 60 sec) , while the large fluctuations of the middle trace are due to the superposition of 
a fluctuating number of events (r = 2600 sec). For the bottom trace, r = 10000 sec 

values reported for x are much smaller than those for Xi- Moreover, correction meth- 



ods (Duri et al., 2005 ) to remove the contribution of the statistical noise to x are often 
used: using these corrections, one has x — for homogeneous dynamics. 

While time resolved correlation (TRC) experiments probe temporal fluctuations of 
the dynamics, they still lack spatial resolution. In photon correlation imaging (PCI), 
by contrast, spatial resolution is achieved by modifying the collection optics. As shown 
in ( |Duri et al., 20096 ) and sketched in Fig. 10.31 for 9 = 90°, a lens is used to image the 
scattering volume onto a CCD or CMOS detector, while a diaphragm limits the range 
of q vectors accepted by the detector. Under these conditions, each pixel of the sensor 
is illuminated by light issued from a small region of the sample and scattered in a 
small solid angle associated with the same, well defined scattering vector. The images 
are analyzed in the same way as for TRC, except that they are divided in regions of 
interest (ROIs): a local degree of correlation ci(r,t,T) is calculated for each ROI, by 
averaging the intensity correlation function over a small set of pixels centered around 



xviii 



■ ■ : ■ 


» - 














.,.=,. 


■ 






• - 


I 


E 










- 




- 




> 
























■ 












■ 








* 

r 


} E* 




- 












- 


_ 


■ 1 - _ 




■ 




















1 


_ 


















- 






; ■ & 


■r 






■ - .- 


■ j' 




. t. - - 








— — 

■ * k-s 




" 
















... * T"' 




■ : 












■■- Jr . '"- 




"V; 


■ 








— 


1 ^ 






■ 










r 

> 


■ 












■ 




- 


- - "- . 


' - 1 







Fig. 0.3 left: scheme of the Photon Correlation Imaging apparatus for a scattering angle 
9 — 90° . The lens L makes an image of the sample S onto the CCD detector. The diaphragm 
D, placed in the focal plane of L, selects only light rays scattered at 9 m 90°. Right: typical 
CCD image recorded by a PCI apparatus. The overlaid boxes indicate the grid of ROIs for 
which local degrees of correlation, ci(r,t,r), are calculated. The size L of the ROIs has been 
exaggerated for the sake of clarity; typically, L is of the order of 10-20 pixels, corresponding 
to ~ 20 — 100 ^im in the sample 

r. The spatial correlation of the dynamics can then be measured by comparing the 
temporal evolution of cj(r, t, r) for ROIs separated by a distance Ar. More specifically, 



we define (Duri et al, 200961 



94 (Av) = B(r)( - <^(M,T)fe j( r + Ar,t,T)) t \ 
\ \J var [<5c/(r, i, r)]var[<5c/(r + Ar, t, t)] j 

where Sci = cj— < Ci >t are the temporal fluctuations of the local dynamics and 
B(t) is a normalizing coefficient chosen so that ^(Ar, r) — > 1 as Ar — > 0. This is the 
analogous, albeit at a coarse grained level, of the spatial correlation of the dynamics 
calculated in numerical and experimental work where particle trajectories are accessi- 
ble (see Chapters 2 and 5). In most cases, the dynamics are isotropic and is averaged 
over all orientations of Ar. It is important to distinguish the length scale over which 
the dynamics are probed from the spatial resolution with which the local dynamics can 
be measured. The former is dictated by the inverse scattering vector. Depending on the 
scattering angle, typical values range from a fraction of up to tens of /Um. The lat- 
ter is determined by the size of the ROIs and the magnification with which the sample 
is imaged. Typical values are in the range 20-100 fim. As a final remark, we note that 



the differential dynamic microscopy method (Cerbino and Trappe, 2008) mentioned 
at the end of Sec. 10.2.11 could be easily adapted to calculate g±. This would improve 
the spatial resolution as compared to that of PCI, thanks to the larger magnification 
typically used in a microscope. 
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Fig. 0.4 Intermediate scattering functions for colloidal hard spheres ([Brambilla e t al, 2009| 
|Brambilla et al, 2010[ ). Data are labeled by the volume fraction ip. The lines are stretched 
exponential fits to the final decay of fs(q,r). Adapted from (Brambi lla et al, 2009[ ) with 
permission 

0.3 Average dynamics and dynamical heterogeneity in the 
supercooled regime 

0.3.1 Structural relaxation time 

The average dynamics of colloidal hard spheres in the supercooled regime (ip < tp g ) has 
been thoroughly studied in a series of works on PMMA-based systems (Pusey and van Megen, 1987 



|van Megen and Underwood, 1994[|"van Megen et al., 1998| [Brambilla et al, 2009[|E1 Masri et al, 20090 . 
Figure 10.41 shows typical ISFs measured for a variety of volume fractions at a scat- 
tering vector q = 2.5/et (a = 100 nm) ( |Brambilla et al, 2009} |Brambilla et al, 20T0| , 
below the first peak of the static structure factor [we quote here the more precise 



determination of a reported in (Brambilla et al, 2010), slightly smaller than that 



in QBrambilla et al, 2009[ )1 . These experiments are performed close to the best index 

matching conditions for a PMMA sample with size polydispcrsity a = 12.2% (Brambill a et al, 2010[ ). 

Under these conditions the sample is optically polydisperse, as discussed in Sec. 10.2.31 

thus, the self part of the ISF, f s , is probed ( |E1 Masri et al, 2009 ). At low volume 

fractions, the decay of f s is well fitted by a single exponential, as expected for diluted 

Brownian particles. As tp increases, the ISFs develop a two-step relaxation. The initial 

decay depends weakly on tp and corresponds to the motion of a particle in the cage 

formed by its neighbors. The final decay corresponds to the relaxation of the cage; its 

characteristic time, r Q , increases by almost 7 decades in the range of (p investigated, 

where all samples equilibrate. 

Figure 10751 shows r Q (ip), as obtained by fitting the final relaxation of the ISF to a 
stretched exponential: 

f s = Bcxp[-(r/T a f}, (0.9) 

with /3 w 0.56 in the glassy regime. In the range 0.517 < (f < 0.585, corresponding to 
about three decades in relaxation time, the volume fraction dependence of r Q agrees 
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Fig. 0.5 Main panel: comparison of the volume fraction dependence of the structural re- 
laxation time, self diffusion coefficient and low shear viscosity as determined in various 
works on PMMA colloidal hard spheres with size polidispersity ranging from » 4% to 
12.2%. All quantities are normalized with respect to their value in the tp — > limit. Solid 
squares fBrambilla et al., 2009[ l, solid circles ( |van M cgcn and Underwood, 1994) and solid 
diamonds ( |van Megen et al., 1998[ l are T a data obtained by DLS. Semifilled symbols in- 
dicate data where only a partial decay of f(q, r) could be measured. Stars are self dif- 
fusion data obtained by DLS ( |Segre et al., 1995c[ ). Open squares and open triangles are 
viscosity data from Refs. ( |Segre et al, 1995 c| and ( |Cheng et al, 2002[ ), respectively. Inset: 
same data as a function of scaled volume fraction, where scaling factors ranging from 1 
to 1.05 were chosen so as to superimpose all curves for tp < 0.2. [The data sets from 
Refs. ( |van Megen and Un derwoo d, 1994] |van Megen et al, 1998[ ), for which no \ow-tp data 
are available, are not included in this plot]. A reasonably good collapse is obtained also for 
ip > 0.2, suggesting that the main source of discrepancy between the various data sets lies in 
the uncertainty on the absolute volume fraction determination 

very well with the critical law predicted by MCT flGotze, 1999[ ): 




(0.10) 



with 7 = 2.6 and (p c = 0.59. This is shown in Fig. 10.6a . where r Q is plotted against 
(1 — ip/ipc)^ 1 . Deviations are observed at low tp (as expected, since here MCT does 
not apply) and, more importantly, at the highest volume fractions that could be 
probed, where r Q grows much slower than expected from MCT. This suggests that 
the divergence predicted by MCT is in fact avoided, as confirmed unambiguously 
by the fact that equilibrium ISFs could be measured up to tp = 0.598, above the 
critical packing fraction tp c = 0.59 obtained from the MCT fit (Brambi lla et al., 2009| 
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Fig. 0.6 a): double logarithmic plot of the structural relaxation time r a as a function of 
(1 — if/ifc)^ 1 , where the critical volume fraction ip c is obtained from a fit to Eq. (|0.10[) 
(adapted from (Brambill a et al, 2009[ ) with permission). In this representation, the MCT law 
is a straight line with slope 7 = 2.5 (solid line). Data at tp > ip c are not represented in this 
plot, b): data from Ref. (Brambil la et al, 2009[ ) for tp > 0.41, plotted using reduced variable 
such that the generalized VFT law, Eq. (|0.11[) with <po — 0.637 and 5 = 2, corresponds to the 
straight line shown in the plot (reproduced from (Brambill a et al, 2009[ ) with permission) 



|E1 Masri et al, 2009[ ). Fi gurc l0.5l shows also data taken from previous works on PMMA 
systems ( |van Megen and Underwood, 1994| |van Megen et al, 1998[ ) , which exhibit a 
similar (apparent) algebraic divergence of r a . Note however that in these works data at 
very high ip could not be obtained, possibly explaining why the crossover between the 
MCT regime and the high tp "activated" regime reported in (Brambil la et al, 2009| 
|E1 Masri et al, 2009[ ) was not observed previously. 

A similar crossover is observed in molecular glass formers ( |Donth, 2001| ): the re- 
duced dynamical range accessible in colloids is responsible for the difficulty in pinpoint- 
ing it, since typically ip c is close to tp g , the volume fraction above which equilibrium 
dynamics become too slow to be experimentally accessible. The analogy with molec- 
ular glass formers suggests that T a (tp) may be fitted by a Vogel-Fulcher-Tammann 
(VFT)-like form: 



, cxp 



.4 



(0.11) 



where the exponent S has been introduced for the sake of generality (6=1 for the 



VFT law). Although a very good fit of the data of Ref. (Brambilla et al, 2009) is 
obtained for S = 1, a somehow better fit is obtained for 5 = 2 and ipo = 0.637. The 
quality of the fit thus obtained can be appreciated in fig. 10.6b . where the data of 
Ref. ( jBrambilla et al, 2009[ ) are plotted using reduced variables so that Eq. (|0.11[) 
reduces to a straight line. A crucial question is whether <pg should be identified with 
<PRCP- Experimentally, this question is still open, due to the difficulties in measur- 
ing precisely the volume fraction at random close packing. Numerical simulations for 
a binary mixture of hard spheres (Berthier and Witten, 20091, however, show that 
ipo = 0.641 < <y9R,cp < 0.664 supporting the thermodynamic glass transition scenario 
and implying that in colloidal hard spheres the glass and the jamming transition are 
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distinct phenomena flKrzakala and Kurchan, 2007 Berthier and Witten, 2009 1 . This 
viewpoint is still highly debated (see, e.g., ( Xu et al., 2009 Kami en and Liu, 2007D ). 

As a final remark on the structural relaxation in colloidal hard spheres, we note that 
one may wonder whether size polydispersity may have a significant impact on the dy- 



namical behavior (van Megen and Williams, 2010 Brambilla et al., 20101 and in par- 
ticularly on the location of the (apparent) MCT divergence. Although no systematic 



experiments have been performed to address this issue ( Williams and van Megen, 2001 ) , 
computer simulations ( |E1 Masri et al., 200 9 Pusc y et al., 20 09; Zaccarc lli et al., 2009a[ ) 
have shown that <p c is essentially unaffected by a in the range 3% < a < 12% typically 
explored in experiments. 



0.3.2 Viscosity 

Viscosity measurements arc an alternative way to probe the average dynamics of hard 
sphere suspensions. Several works have been devoted to the <p dependence of the viscos- 



ity, t\, in the limit of a vanishingly small applied shear, see e.g. Ref. ( Cheng et al., 2002 1 
and references therein. These measurements are quite delicate, since the applied stress 
has to be very small (as low as a fraction of a mPa) in order to avoid non-linear ef- 
fects, and because the resulting shear rate becomes extremely small beyond ip w 0.5, 
limiting measurements in the deeply supercooled regime. Additionally, the comparison 
between results obtained for different samples is affected by the uncertainties on the 
absolute volume fraction discussed above. Cheng et al. ( |Cheng et al., 2002 ) show that 
data from various groups collapse reasonably well on a master curve when the absolute 
volume fractions are scaled by a factor up to about 1.03 to account for polydispersity. 
Figure [Q~5l shows viscosity data from Rcfs. (C heng et al., 20 02; Scg re et al., 1995c| to- 
gether with the DLS data discussed in Sec. 10.3.11 While there is some discrepancy 
between DLS and viscosity data, these differences are likely to be due, to a great 
extent, to uncertainties in ip and to the different methods used in the determination 
of the absolute volume fraction. Indeed, data sets for which data points at low ip are 
available can be scaled reasonably well onto a master curve by correcting ip using 
scaling factors between 1 and 1.05, so as to superimpose the growth of the viscosity 
or the relaxation time at low volume fraction (ip < 0.2) (see the inset of Fig. |075|) . The 



relationship between r/ and r a is discussed in more detail in Ref. (Scgre et al., 1995c 



where it is found that the low shear rate viscosity and the structural relaxation time 
measured by the collective ISF at the peak of the structure factor agree remarkably 
well up to ip ~ 0.5, while the relaxation time for the self part of the ISF is somehow 
lower. 

The nature of the divergence of rj is still a matter of debate. In their early work ( |Phan et al., 1996 ) 
Russel, Chaikin and coworkers reported that r/((p) is fitted well by the Krieger-Dougherty 
equation, a critical law of the same form as Eq. (|0.10[) with 7 = 2, yielding ip c = 0.577, 
consistent with the critical packing fraction obtained from MCT fits of DLS data. How- 
ever, in their subsequent analysis of data on a larger range of viscosity ( |Cheng et al., 2002] ) , 
they report that a VFT-like fit (e.g. of the form of Eq. (|0.11|) with 6 = 1) yields better 
results. They find po = 0.625, close to random close packing, and significantly higher 
than tp c . It should however be recalled that viscosity measurements are feasible only 
up to volume fractions lower than <p c (e.g. ip < 0.562 in Ref. QCheng et al., 2002] )), 
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making it difficult to draw unambiguous conclusions on the nature of the divergence 
of ?], and in particular on the existence of a divergence at <p c « 0.58. 

0.3.3 Dynamical heterogeneity 

The first microscopy experiment that examined the motion of supercooled colloidal 
particles was by Kasper, Bartsch, and Sillescu in 1998 ( [Kasper et al, 1998[ ). They 
devised a clever colloidal system primarily composed of refraction-index-matched par- 
ticles made from cross-linked poly-t-butylacrylate. They then added a small concen- 
tration of tracer particles which had non-index-matched cores of polystyrene coated 
with shells of poly-i-butylacrylate. Using dark field microscopy, they could observe 
the motion of the tracer particles. They observed that particles exhibited caged mo- 
tion, as described above in Sec. 10.3.11 That is, a particle would diffuse within some 
local region, trapped in a cage formed by its neighbors, and then occasionally ex- 
hibit a quicker motion to a new region. This was useful evidence that the dynamics 
are temporally heterogeneous, and the first direct experimental visualization of caged 
motion. Averaging over all of the tracer particles, they noted that the distribution 
of displacements was non-Gaussian, likely linked to the cage trapping and cage rear- 
rangements. This experiment exploited an inherent size polydispersity of about 8% to 
prevent crystallization. The chief limitations of the experiment were that the observa- 
tions were limited to two-dimensional slices of the three-dimensional sample, and also 
that only isolated tracer particles were observed, rather than every particle. 

The next published experiment, by Marcus, Schoficld, and Rice, had different trade- 
offs flMarcus et al, 1999[ ). They used a very thin sample chamber to study a quasi- 
two-dimensional colloidal suspension; the spacing between the walls of their sample 
chamber was approximately 1.2 particle diameters. Because of the relative ease of 
studying a thin sample, they did not need tracer particles, but rather could follow the 
motion of every particle within the field of view. Like Ref. ( |Kasper et al, 1998D , they 
observed cage trapping and cage rearrangements, and found a non-Gaussian distribu- 
tion of displacements. Being able to see all of the particles, they also noted that the 
cage rearrangement motions were spatially heterogeneous, with groups of particles ex- 
hibiting string-like motions. The string-like motions were quite similar to those seen in 
simulations < \Kob et al, 1997[|Donati et al, 19981 |Hurley and Harrowell, 19960 . Their 
results clearly showed a connection between the non-Gaussian behavior and the spa- 
tially heterogeneous dynamics, in that the non-Gaussian displacements were due to 
the particles involved in the cage rearrangements. 

The experiment of Kegel and van Blaadcrcn in 2000 used confocal microscopy to im- 
prove upon the prior experimental limitations flKegel and van B laadcrc n, 2000] ). They 
observed the motion of core-shell colloidal particles, in a fully three-dimensional sample 
(although their observations were limited to two-dimensional images to maximize the 
imaging rate). Their colloidal samples were well-characterized as having hard-sphere 
interactions. In this experiment, they again observed string-like regions of high mobil- 
ity, related to the non-Gaussian distribution of displacements. This was the first ex- 
periment to directly visualize spatially heterogeneous dynamics in a three-dimensional 
sample, confirming what simulations had already suggested, that string-like motion 
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Fig. 0.7 Picture depicting the current positions of the most mobile colloidal particles (light 
colored) and the direction they are moving in (dark colored) . The particles are drawn at 75% 
of their correct size, and for clarity only the most mobile particles are shown. Many particles 
move in similar directions to their neighbors, as shown by those within the oval. However, 
some particles also move in opposite directions to their neighbors, for example closing gaps 
between them, as highlighted by the rectangle. This is a sample with ip — 0.56, data taken 
from Ref. ( | Weeks et al, 2000[ ). The particles have radius a — 1.18 fj,m and the time scale 
used to determine the mobility is At = 1000 s 



is not an artifact of two dimensional systems (Marcus et at, 1999 Kob et al., 1997 



Donati et al., 1998 Hurley and Harrowell, 1996 



Perera and Harrowell, 1999). 



Shortly after Ref. ( Kegel and van Blaaderen, 2000 ), Weeks et al. published a sim- 
ilar experiment using confocal microscopy to study three-dimensional colloidal sam- 
ples ( |Weeks et al., 20 00). Utilizing a faster confocal microscope than Ref. ( [Kegel and van Blaaderen, 2000[ ), 
they were able to observe displacements in three dimensions. They too saw string-like 
motion, although also noted that some particles were moving in non-string-likc ways 
termed "mixing " d Weeks and Weitz, 2002] ); see Fig. [0/7] A limitation of this experi- 
ment is that the particles were later discovered to be slightly charged, rather than being 
ideal hard spheres ( |Gasser et al., 200T| ). The three-dimensional observations enabled 
the fractal nature of the regions of mobile particles to be measured as df = 1.9 ± 0.4, 
similar to simulations (Donati et al., 19991. The sizes of the mobile regions increased 



dramatically as the colloidal glass transition was approached. 

Subsequently, the data of Ref. ( | Weeks et al., 2000] ) has been reanalyzed to highlight 
other features. The idea of caging was quantified in Ref. ( Weeks and Weitz, 2002[ ) , find- 
ing a result similar to Kasper et al. ( |Kasper et at, 19"98 ), that caging manifested itself 



as an anti-correlation of particle displacements in time. That is, if a particle pushed 
against the "walls" of its cage (formed by its neighbors), then the particle was likely 
to subsequently be pushed backwards. The original work of Ref. ( | Weeks et al., 2000] ) 
highlighted mobile regions using a slightly arbitrary definition of which particles were 
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mobile. Later analysis used spatial correlation functions to avoid defining a particular 
subset of particles as mobile ( |Weeks et ai, 20 07). These correlation functions found 
length scales for particle mobility which grew by a factor of 2 as the glass transition was 
approached, finding the largest length scale for a super-cooled liquid of approximately 
8a in terms of the particle radius a. 

More recent work has taken the study of dynamical heterogeneities in colloidal 
suspensions in new directions. One notable set of experiments uses superparamagnetic 
particles and an external magnetic field to control the glassiness of a sample in situ 



( |Konig et al., 2005||Ebert et ai, 2009||Mazoyer et ai, 2009 1. Other experiments study 
the colloidal glass transition in confinement ( |Nugent et ai, 2007 ), or glassy samples 



as they are sheared ( |Besseling et al, 2007] |Schall et al, 2007[ |Chen et al., 2010D 



Dynamical heterogeneity in the equilibrium regime has also been probed by DLS. 
Attempts to directly measure x± or 94 using the time- and space-resolved methods 
discussed in Sec. 10.2.41 have, so far, failed, because these techniques lack the reso- 
lution needed to detect dynamical heterogeneity on the length scale of a few par- 
ticles. By contrast, Xi can be measured indirectly using the theory developed in 
Chapter 3 (see Sec. 3.2.5 therein). For colloidal hard spheres the following relation 
holds ( [Berthier et ai, 2005| |Berthier et al, 2007fe[ |Berthier et al., 2007a[ ) : 

X4,(q,r) = Xi(q,T)\ v + pk B TK T [ipx v (q,T)} 2 , (0.12) 

where p is the number density, kt the isothermal compressibility (taken from the 
Carnahan-Starling equation of state), Xii.Qj T )\v denotes the value taken by X4(<3S T ) m 
a system where density is strictly fixed, and x^( ( ?) T ) = df s (q, r)/dip. Only the second 
term in the r.h.s. of Eq. (|0.12[) can be accessed experimentally. Numerical simulations, 
where both terms in the r.h.s. of (|0. 12[) can be calculated ( |Brambilla et ai, 2009[ ), 
show that the first term can be neglected in the deep supercooled regime: X4(q, r) w 
pkBTnTlcfXipili T )} 2 when the latter term is larger than unity. Experimentally, x<p 
can be obtained either by numerical differentiation using two ISFs measured at close 
enough volume fractions (see Fig. 10.8a ). or by using the chain rule in the r.h.s. of 
Eq. dnU): 

df, = df 1 dB + df ± dr^ + df 1 dl 
dtp OB dtp dr a dtp d(3 dip' 

where the partial derivatives with respect to volume fraction of the coefficients B, 
r a and (3 defined in Eq. (|0.9p are obtained by fitting their ip dependence by smooth 
polynomials ( |Dalle Ferrier et al., 2007] |Brambilla et al., 2009[ ) . 

Figure ID~5b shows X4{q = 2.5a, r) together with the ISFs used to obtain it by finite 
difference ( |Berthier et ai, 2005 ), for the same system as in Ref. (Brambilla et al., 20091 



The dynamical susceptibility has the characteristic peaked shape observed in a variety 
of glassy systems (see also Chapters 2 and 5), with the maximum of the fluctuations, 
X* = Xi( T *)^ occurring on a time scale r* comparable to the structural relaxation 
time r Q . As discussed in detail in Chapters 2 and 5, x* is of order iV corr , the number of 
particles that undergo correlated rearrangements. Figure [Ch8b shows x*, obtained for 
the same system but using the chain rule (|0.13j) . as a function of <p. The amplitude of 
dynamical fluctuations increases with volume fraction, supporting theoretical scenarios 
where a growing dynamical length accompanies the divergence of the relaxation time 
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Fig. 0.8 a), left axis: ISFs at two nearby volume fractions (ip = 0.5953 and 0.5970 for the 
open squares and the crosses, respectively). The solid circles are \ v as obtained by finite 
difference from the two ISFs (right axis, same scale as in b)). b): height of the peak of X4> 
X* , as a function of ip (adapted from ([Brambill a et al., 2009 [ ) with permission). Here, \a has 
been estimated using the chain rule, as explained in the text. The dashed line is a MCT fit to 
the dynamical susceptibility. The vertical dotted line indicates the location of the divergence 
predicted by MCT but avoided by the data 

on approaching a glass transition. The dashed line shows the predictions of advanced 
mode coupling theories, x* ~ (<Pc — tp)~ 2 (Bcrthic r et al., 2007fc| Berthier et al., 2007a 
Biroli and Bouchaud, 2004~| . As observed for the average dynamics, the data agree with 
MCT only over a limited range of volume fractions; in particular, MCT overprcdicts 
the growth of \* & t high ip, where the data remain finite across the the critical volume 
fraction ip c . Therefore, the analysis of dynamical fluctuations confirms the crossover 
from an MCT regime to an activated regime at high volume fractions inferred from 
the average dynamics. 



0.4 Average dynamics and dynamical heterogeneity in 
non-equilibrium regimes 

0.4.1 The glassy regime 

Thus far we have focused on dynamical heterogeneities in supercooled colloidal liquids; 
these systems will crystallize after some time (if the particles are sufficiently monodis- 
perse). For polydisperse systems, the dynamics can be stationary (unless the sample 
has been recently sheared or stirred), since the sample is in (metastable) equilibrium. 
In contrast, glasses (volume fraction ip > ip g ) are out of equilibrium, and their prop- 
erties depend on time, a phenomenon termed "aging." In particular, consider a low 
volume fraction colloidal suspension which is centrifuged to rapidly increase its volume 
fraction to the point where it becomes a colloidal glass at some time t w = 0. The mo- 
tion of particles within this sample depends on the time t w since this formation, called 
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the waiting time or simply the age of the sample. Initially, particles can move rela- 
tively rapidly, but as the system ages, particle motion slows ( |van Megen et ai, 1998| ) . 
Particles take longer to move the same distance that was covered quickly at an earlier 
age. Equivalcntly, the cages formed by a particle's neighbors are more long-lasting, 
and given that these samples have a high volume fraction, particles spend almost all 
of their time tightly confined within their cages. 

In the earliest microscopy experiments (Ref. ( |Kasper et al, 1998[ )), there was evi- 
dence of slow motion within glassy samples. While the mean square displacement of 
the particles grew extremely slowly, it did grow. This provided some evidence that 
particles were not completely frozen, but rather moved to new positions as the sample 
aged. 

More direct observations were obtained via confocal microscopy in Ref. ( | Weeks et ai, 2000| . 
In that experiment, Weeks et al. observed that while mobile regions of particles grew 
larger as tp — > (p g (in the super-cooled state), in glassy samples the mobile regions 
were quite small. This implied that spatial dynamical heterogeneities were less impor- 
tant for colloidal glasses. This result was revised in 2003 by Courtland and Weeks, 
who observed larger clusters of mobile particles in an aging colloidal glass, again us- 
ing confocal microscopy ( |Courtland and Weeks, 200 31- The key difference was in the 
data analysis. In a glassy sample, most of the particle motion is Brownian motion 
within the cages. Occasionally particles have cage rearrangements, but this is hard to 
distinguish from the Brownian motion as the distances particles move during these re- 
arrangements gets very small fl Weeks and Weitz, 2002| |Courtland and Weeks, 200"3| . 
Courtland and Weeks were able to observe the cage rearrangements by low-pass fil- 
tering the raw trajectories, thus smoothing out the Brownian motion and making the 
cage rearrangements clearer. They found that there was indeed spatially heterogeneous 
dynamics comparable to the behavior of super-cooled colloids, but surprisingly they 
did not see any dependence on the aging time t w . Similar results were also seen in a 
later confocal microscopy study of a binary colloidal glass ( |Lynch et ai, 2008| . 

One important consideration for studies of the aging of colloidal glasses is the 
protocol for achieving the initial glassy state. In traditional molecular glasses, a sam- 
ple is quenched from a liquid state by rapidly reducing the temperature. In a col- 
loidal glass, the analogous quench would be to rapidly increase the volume frac- 
tion, for example by centrifugation. However, in the experiments described above 
( |Courtland and Weeks, 2003| |Lynch et ai, 20081 ) i the samples were instead prepared 
at a constant volume fraction, and then shear-melted by stirring them with an em- 
bedded stir bar. This is sometimes termed "shear-rejuvenation" as the stirring will 
make a well-aged colloidal glass look like a young colloidal glass, that is, it resets 
t w = 0. However, there is evidence that these two protocols give different glassy states 
in molecular glassformers (McKcmi a, 2003p . 

Motivated by this, two groups have found ways to quench a colloidal glass in situ 
from a low volume fraction state to a high volume fraction state. One method uses 
an external magnetic field to control the effective inter-particle attraction in a two- 
dimensional colloidal sample composed of superparamagentic particles ( |Assoud et al., 2009[ ). 
This technique has not been used to study aging explicitly, but so far has focused on 
slow crystallization after a rapid quench ( |Assoud et ai, 200 9). Another method uses 
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the sample temperature to control swelling in hydrogel particles, again in a quasi- 
two-dimensional experiment ( Yunker et ai, 2009[ ). Experiments using soft swellablc 
particles are described in Sec. 



(TO 



Dynamical heterogeneity in glassy colloidal samples have also been probed by 
DWS. We recall that this light scattering technique is sensitive to motion on very 
small length scales, down to a fraction of a nm (see Sec. 10 . 2.3)) . a highly desirable fea- 
ture for glassy systems where particles hardly move. Reference QBallesta et ai, 2008 § ) 
discusses both the average dynamics and its temporal fluctuations in very dense sus- 
pensions of relatively large particles (a f=s 10 /mi). After initializing the sample by 
shaking it vigorously, the dynamics slow down until a pseudo-stationary regime is at- 
tained, where all measurements are performed. This regime does not correspond to a 
true equilibrium state, but rather to a regime where the very local dynamics probed 
in this experiment do not evolve significantly on the time scale of the experiments (up 
to a few days). Figure IUT91 shows the volume fraction dependence of various parame- 
ters characterizing the dynamics in the pseudo-stationary regime. As shown in a), the 
average relaxation time of the intensity correlation function, as obtained from a fit 
92{t,r) — 1 = £?exp[— (t/to^))^*)], increases smoothly with ip, seemingly diverging 
at ip = ip max = 0.752, presumably the random close packing volume fraction for this 
quite polydisperse system. Panel b) shows the ip dependence of the average stretching 
exponent, j3, which increases from about 0.9 to 1.3 as random close packing is ap- 
proached. Panel c) shows the height, x* of the peak of the dynamic susceptibility x 
defined in Eq. (|0.7[) . Quite surprisingly, x* is non-monotonic: it first increases with ip, 
reaches a maximum, but eventually drops dramatically close to random close packing. 

The interpretation of these data requires special care, since the usual DWS formal- 
ism has been developed for spatially and temporally homogeneous dynamics. By con- 
trast, Ref. flBallesta et ai, 2008 fc| ) discusses a general model for DWS for heterogeneous 
dynamics, based on an extension of the formalism originally proposed by Durian et al. 
for DWS measurements of the spatially localized, temporally intermittent dynamics 
of a foam ( Durian et ai, 1991[ ). According to the model of (Ballesta et al., 200861, the 
non-monotonic behavior of x* (y) results from the competition between two contrast- 
ing effects. On the one hand, the size £ of the regions that undergo a rearrangement 
increases continuously with volume fraction, leading to enhanced temporal fluctua- 
tions as observed in supercooled samples (see Sec. 10 . 3.3|) . On the other hand, as <p 
increases the particle displacement associated with each of these events is increas- 
ingly restrained, due to crowding. Thus, an increasingly large number of rearrange- 
ment events is required to decorrelate the scattered light at high ip, leading to re- 
duced fluctuations. While simple simulations of a DWS experiment in a medium 
undergoing rearrangements described by this scenario reproduce the experimental 
data (Ballest a et ai, 2008fc[ ), a direct measurement of £(</?) would be necessary to con- 
firm it. This would allow one to better understand analogies and differences with the 
experiments on driven grains (Lechenault et al, 2008) mentioned in Chapter 5, where 
a non-monotonic behavior of both x* and £ has been reported. 
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Fig. 0.9 Volume fraction dependence of the dynamics and its temporal fluctuations for 
a concentrated suspension of colloids as probed by DWS. Solid symbols refer to fresh 
samples, semi filled symbols to aged samples that have been rejuvenated mechanically. 
See ( |Ballesta et al, 2008b| ) for more details, a) Average relaxation time, the line is a crit- 
ical law fit to the data, tq ~ (y>max — <p)~ y with </? m ax = 0.752 and y — 1.5. Here and in the 
other panels the vertical dashed line indicates the location of fi mla . b) Stretching exponent 
obtained by fitting the final relaxation of the intensity correlation function to a stretched 
exponential, c) Height \* of the peak of the dynamical susceptibility, \(t). Note the abrupt 
drop of x* on approaching (fimux- In a) and b), the bars indicate the standard deviation of 
the temporal distribution of the relaxation time and the stretching exponent, respectively. 
Adapted from Ref. ( |Ballesta et al, 20086| > with permission 

0.4.2 Dynamical heterogeneity under shear 



Aging systems are out of equilibrium; another type of non-equilibrium system is a 
driven system. Of particular interest are samples that are sheared. The importance 
of shear can be quantified by a Peclet number (see also Sec. I0.1.2[) . This is the ratio 
of the time scale for diffusion to the time scale for shear-induced motion. The shear 
induced time scale is given in terms of the strain rate as 1 /j. Recall that the Brownian 
time scale is tb = a 2 /(6D) (Eq. (j0.1|l ). Typically when considering dense colloidal 
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suspensions, the relevant diffusion constant is D 00l the long-time diffusion constant, 
which varies with the volume fraction. The other option would be to use Dq, as given 
by the Stokes- Einstein-Sutherland formula, Eq. (|0.2[) . which is the diffusion constant 
in the (p — > limit. If Do is used, the Peclet number is termed the bare Peclet number, 
Pe, and if is used, it is termed the modified Peclet number, Pe* . Given that 
we wish to understand how particles move and rearrange, this is the long-time-scale 
motion (Doo), and combining the expressions above we find Pe* = a 2A //(6D oc ). For 
Pe* < 1, diffusion is the primary influence on particle motion, and the shear is only 
a small perturbation. For Pe* > 1, the shear- induced motion is expected to be more 
significant. The main consideration here is that for the shear rate to be significant, an 
experiment must have Pe* > 1, which is the case for the experiments described below 
in this subsection. 

Simulations of supercooled fluids found that as Pe* increases, the system becomes 



more-liquid like, and dynamically heterogeneous regions become smaller ( Yamamoto and Onuki, 1997), 
similar to the idea of shear unjamming a sample ( |Liu and Nagel, 1998[ ). As yet, this 
relationship between Pe* and dynamical heterogeneity has been untested by colloidal 
experiments. 

A complementary question is to ask what the nature of a shear-induced motion 
is, for a dense amorphous sample. For example, sheared crystalline materials respond 
by the internal motion of dislocation lines. Influential simulations by Falk and Langer 
found, for amorphous materials, that rearrangements occurred in localized regions, 
termed "shear transformation zones" ( |Falk and Langer, 1998[ ). In these regions, locally 
the stress builds up until it is released by a rapid local rearrangement event. These 
behaviors have been directly observed in an experiment by Schall et al., where they 
used confocal microscopy to examine a colloidal glass as it was sheared between two 
parallel plates at low strains ~ 4% flSchall et al., 200 7). In their experiment, the shear- 
induced dynamical heterogeneities had a small spatial extent, although in some cases 
they appeared to relax the strain over a large region even though the farther away 
particles did not move very far. Their activation energy was estimated to be E ~ 
16fcsT, where fceT is the thermal energy. Given the low strain, it is possible that 
these observations are of shear-induced aging effects, rather than shear flow. 

A separate experiment examined the shear-induced motion of colloidal particles in 
supercooled colloidal liquids ( |Chen et al., 2010] ). In these less dense samples, localized 
rearranging regions were observed at high strain values. These samples were sheared 
between two parallel plates that moved back and forth with a triangle wave displace- 
ment curve. Strikingly, the shapes of these rearranging regions were isotropic on aver- 
age, and showed no distinction between the velocity direction, velocity gradient direc- 
tion, and the vorticity direction mutually perpendicular to the first two. Similar obser- 
vations were made of a sheared colloidal glass ( |Besseling et al., 2007[ ) : in those experi- 
ments, little difference was found in the effective diffusivity in the three directions, once 
the affine motion of the average shear profile was subtracted from the particle displace- 
ments. This study was the first direction observation of a sheared colloidal glass which 
directly imaged shear- induced dynamical heterogeneities ([Bcsscling et al., 2007[ ) . Note 
that the isotropic nature of plastic rearrangements observed in experiments is at odds 
with recent simulations of a 2D supercooled fluid of soft particles, where anisotropic 
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structural rearrangements were observed ( |Furukawa et al., 2009| ). These observations 
depended on careful analysis methods, which have not yet been applied to experimen- 
tal data from colloids. 

While not directly the same idea as shear-induced dynamical heterogeneity, it is im- 
portant to note that many soft glassy systems - such as colloids - exhibit shear-banding 
flDhont, 1999l|Dhont et al, 2003| [Fielding, 2007||Dhont and Briels, 2008||Fielding et al, 20090 . 
That is, when a large sample is sheared, sometimes the strain is localized near one of the 
boundaries. Within the "shear band" the sample is straining a significant amount and 
plastically deforming, while outside the shear band, the sample is nearly unstrained. 
As the stress must be continuous throughout the sample, this suggests that the ma- 
terial is acting as if it has two states, a low viscosity state in the shear band where 
the sample flows, and an elastic state outside the shear band without flow. Shear 
bands have been noted in colloidal suspensions ( Vermant, 2001| |Chen et al, 2010 



|Besseling et al, 2007[|Ballesta et al., 2008a||Dhont and Briels, 20080, colloidal ge ls ( |Moller et al, 20080 , 
worm- like micelles (|Berret et al, 1997 Olmsted, 1999 Cates and Fielding, 20060, foams 



( Debregeas et al., 2001 Lauridsen et al., 2004 Janiaud et al, 20060 , and granular ma- 
terials ( Losert et al., 200'o[ Utter and Behringer, 20080 . Note that the experiments dis- 
cussed in the previous paragraph (Ref. ( |Besseling et al., 2007||Chen et al., 20100 ) were 
observations of a homogeneously shearing subregion within the shear band. 



0.5 Beyond hard spheres 

0.5.1 Soft particles 

Colloidal particles with soft repulsive potential interactions can be obtained in several 
ways. One possibility is to exploit the softness of an electrostatic or magnetic repulsive 
potential. Another possibility is to modify the synthesis of the particles, e.g. in star 
polymers or microgcl particles, where the degree of softness is governed by the number 
of arms and the degree of crosslinking, respectively. Finally, systems based on the self- 
assembly of amphiphilic molecules can form soft spheres, e.g. with amphiphilic diblock 
copolymers. 

For soft systems, a nominal volume fraction, </J n0 m; is often defined, based on 
the size of isolated particles and their number concentration; since particles can be 
squeezed, f> nom can exceed unity. At low concentration, soft spheres behave similarly 
to hard spheres, provided that an effective size that takes into account the range of the 
soft repulsion is used. This is shown, e.g., by viscosity measurements ( |Roovers, 1994[ 
Buitenhuis and Forster, 1997| |Senff et al., 1999 ), where the increase of the relative 



viscosity with concentration is essentially indistinguishable from that of hard spheres. 
At higher concentration, however, the behavior deviates significantly ( Roovers, 1994| 



|Buitenhuis and Forster, 1997||Sessoms et al., 2009; M attsson et al., 2009 00 : the growth 
of the viscosity or the relaxation time with ip nom is much gentler and samples with 
Vnom > </5rcp may still be fluid, since particle deformations allow for structural relax- 
ation. 

In the past years, aqueous solutions of poly-N-isopropylacrylamidc (PNiPAM) mi- 
crogels are become one of the most popular soft systems. Not only can their soft- 
ness be controlled during the synthesis (by varying the amount of crosslinking), but 
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it can also be tuned by varying the temperature, T. Under appropriate conditions, 
a decrease of a few degrees of T results in a ~ 20% growth of the particle ra- 
dius (Senff et ai, 19991 and in an increased softness. The influence of the softness 
of concentrated PNiPAM particles on their dynamics has been studied in detail by 
DLS in Ref. (Mattsso n et ai, 2009 djj . Quite remarkably, the authors find that soft- 
ness correlates with fragility, defined here as the slope of logTc vs concentration, 
£, at the glass transition, in analogy with the definition for molecular glass form- 
ers ( Donth, 2001[ ), where 1/T replaces £. Very soft particles have an Arrhenius-like 
behavior, r a ~ exp(C£), harder particles have a larger fragility (steeper increase of 
T a(()), and hard spheres are the most fragile system, with the steepest increase of r Q 
on approaching the glass transition. Thus, soft colloids appear as a promising system 
for understanding the origin of fragility in glasses, a long standing problem. 

A similar system has been investigated by photon correlation imaging (PCI) in 
Ref. ( Sessoms et al., 2009 1. Both x and <74 have been studied as a function of (f nom in 
measurements of the collective dynamics at low q. The amplitude of temporal fluctua- 
tions of the dynamics is found to increase monotonically with (p nom , while the range of 
spatial correlations of the dynamics has a non-monotonic behavior, a maximum being 
observed at v'nom ~ ¥>rcp- Remarkably, spatial correlations of the dynamics here ex- 
tend over several mm, in analogy with what was reported for other jammed materials 
(see Sec. lQX2|) . 

The T dependence of the size of PNiPAM particles provides a convenient way to 
quench them rapidly in a glassy state, by preparing a relatively concentrated — yet 
fluid — sample, which is then cooled by a few degrees, thereby swelling the particles 
to a quenched glassy state in a fraction of a second. This protocol was used in a 
quasi- two-dimensional experiment (Yunker et ai, 20091, where observations of aging 
over nearly six decades in aging time t w were possible, due to the rapid quench rate. 
The authors observed spatial dynamical heterogeneity; particles occasionally under- 
went irreversible rearrangements. While the average size of the rearranging regions 
remained approximately constant during aging, similar to the prior study of 3D hard 
sphere samples (Ref. ( |Courtland and Weeks, 2003 1), they identified an increase in the 
size of a particular class of rearranging regions as the sample aged. Specifically, the 
domain size of rearranging particles surrounding irreversible rearrangements increased 
during aging. The largest clusters of rearranging regions involved approximately 100 
particles, a size much smaller than the range of dynamical correlations measured 
in (Sessoms et ai, 20091. However, the technique used to identify rearranging par- 
ticles limited the correlation size, making direct quantitative comparisons difficult. 
Additionally, they saw a relation between the local structure and the particle motion, 
which agrees with some prior work on hard spheres ( |Cianci et al, 2006D - 



0.5.2 Attractive particles 

When colloidal particles experience attractive forces, arrested phases may be obtained 
also at a volume fraction lower than that required for glassy dynamics to be observed 
in hard spheres. The distinction is often made between colloidal gels (up to tp ~ 0.3) 
and attractive glasses (over ip ~ 0.5)( |Trappe and S andkuhlcr, 2004). Very recent work 
by Zaccarelli and Poon further refines the classification of concentrated attractive sys- 
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terns flZaccarelli et al, 2009 fc| ), based on the predominance of either caging or bonding. 
Concentrated, attractive glasses ( |Eckert and Bartsch, 2002"[|Pham et al, 2002[ ) have a 
structure similar to that of repulsive HS systems, while the structure of diluted gels de- 
pends on the strength of the interactions: highly attractive systems tend to form string- 
like, fractal structures, while gel strands are thicker when the interparticlc potential 
well at contact is close to k-gT ([Campbe ll et al, 2005| [Dibble et al, 2006| . The mech- 
anism leading to gel formation is still intensely debated; recent work, still controverted, 
points to the role of an underlying fluid- fluid phase transition which is arrested once the 
dense phase becomes too concentrated ( |Lu et al, 2008[ ) (see Ref. ( |Zaccarelli, 2007[ ) for 
a recent review on colloidal gels). Experimentally, attractive systems are typically re- 
alized either by screening the Coulomb repulsion in charge-stabilized systems, thereby 
exposing the particles to short range, attractive van der Waals forces ( Russel et al., 1992] ), 
or by means of the depletion force ( |Asakura and Oosawa, 1958} |Vrij, 1976] ) induced by 
adding to the suspension smaller particles, often polymer coils ( Poon, 2002[ ). 

Dynamical heterogeneity in weak gels have been explored mainly by confocal mi- 
croscopy ( Gao and Kilfoil, 2007 Dibble et al, 2008 1 and simulations ( |Puertas et al, 2004 



|Charbonneau and Reichman, 2007 ). The general picture emerging from these works 
is that DH is closely related to structural heterogeneity. The probability distribution 
function of the particles displacement (van Hove function) has typically a non-Gaussian 
shape, with "fat" tails corresponding to fast particles whose displacement is anoma- 



lously large (Gao and Kilfoil, 2007 Dibble et al, 2008 Chaudhuri et al, 2008). These 



particles are located at the boundaries of the thick strands constituting the gel, while 
the particles buried within the strands are the least mobile. A detailed analysis of 
particle mobility as a function of the number of their neighbors ( |Dibble et al, 2008| 
confirms this picture. It should be noted that such a structural origin of DH is in 
contrast with hard sphere systems, where no clear connection between DH and struc- 
tural quantities could be established so far. The influence of structure on dynamical 
heterogeneity in attractive systems has also been highlighted in a scries of simulation 
papers by A. Coniglio and coworkers (see e.g. Ref. ( |Fierro et al, 2008"] ) ). 

The length scale dependence of dynamical heterogeneity in weak gels has been 
explored both numerically ( Charbo nneau and Reichman, 200"7| and in XPCS experi- 
ments ( Trappe et al, 2007| . In these works, the peak x* of the dynamical susceptibility 
X4 has been measured as a function of scattering vector q. x* has a non-monotonic 
behavior, the largest dynamical fluctuations being observed on a length scale of the 
order of the range of the attractive interparticlc potential. This has to be contrasted 
to the case of repulsive systems, where the maximum of x* typically occurs around 



the interparticlc distance (Charbonneau and Reichman, 2007 Dauchot et al, 2005). 
On a more technical level, it is worth noting that Ref. (Trappe et al, 2007) has 
demonstrated that modern synchrotron sources and X-ray detectors are now suffi- 
ciently advanced to allow for measurements of dynamical heterogeneities. The activity 

Wandersman et al, 2008[ 



in this field is thus growing rapidly ( Trappe et al, 2007] 



Herzig et al, 2009 Duri et al, 2009a 



Wochner et al, 2009) and there is hope that 



eventually X-ray scattering experiments may probe dynamical heterogeneity in molec- 
ular glass formers and not only for colloidal systems. 

Optical microscopy studies of the dynamics of strong gels are difficult, due to the 
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restrained and very slow motion of particles in these systems flDibble et al, 2008] ) . 
By contrast, scattering techniques have been successfully applied to characterize the 
slow dynamics of tenuous, fractal-like gels made of particles tightly bound by van 
der Waals forces flCipelletti et al, 2000||Duri and Cipelletti, 2006||Duri et al, 2009 60 . 
Time resolved correlation (TRC, see Sec. 10.2 .4[) experiments show that the dynamics 
are due to intermittent rearrangement events where particles move over relatively 
small distances, of the order of a fraction of a fxm flDuri and""C ipcllctt i, 2006] ). Quite 
surprisingly, each of these events affects a macroscopic portion of the sample: the 
spatial correlation of the dynamics measured by photon correlation imaging (PCI, see 
Sec. I0.2.4[) hardly decays over several millimeters, indicating that spatial correlations 
of the dynamics are limited essentially only by the system size ( Duri et al, 20096| . 
Once averaged over both time and space, the intensity correlation function 32(9, t) — 1 
measured by "regular" multispeckle DLS exhibits a peculiar q dependence of both 
the relaxation time and the stretching exponent p obtained by fitting 52 — 1 to a 
stretched exponential, 52(9, t) — 1 ~ cxp[— (t/tv) p ] (Cipelletti et al, 2000). As shown 
in Fig. 10.10a . p is larger than one; accordingly, the relaxation has been termed a 
"compressed" exponential, as opposed to the as opposed to the stretched exponential 
relaxations often observed in glassy systems (p < 1). Moreover, r r ~ q- 1 , as opposed 
to T r ~ q~ 2 as for diffusive motion. 

In Refs. (Cipelletti et al, 2000 Bouchaud and Pitard, 20021 these dynamics were 
interpreted as ultraslow ballistic motion due to the slow evolution of a strain field set by 
internal dipolar stresses. A more refined model has been proposed in Ref. ( |Duri and Cipelletti, 2006[ ), 
taking into account the results from time- and space-resolved light scattering experi- 
ments that highlight the discontinuous nature of the relaxation process. The model is 
based on the following assumptions: i) the dynamics are due to individual rearrange- 
ment events that are random in time (Poissonian statistics); ii) each event affects the 
whole scattering volume (as indicated by PCI); iii) the displacement field induced by 
one single event is that due to the long-range elastic deformation of the gel under 
the action of dipolar stresses (in order to account for p > 1); iv) on the length scales 
probed by the scattering experiments, the displacement due to successive events oc- 
curs along the same direction (i.e. the motion is, on average, ballistic-like, as implied 
by the scaling r r ~ Q 1 )- The model contains just two adjustable parameters: the 
rate of the events in the scattering volume, 7, and the average particle displacement 
resulting from one single event, S. The model captures well the q dependence of both 
p and TV as observed in the average dynamics, as shown in Fig. 10.10a . It also cap- 
tures correctly the growing trend for the q dependency of the amplitude of dynamical 
heterogeneity, although it overestimates their magnitude by about a factor of two, as 
seen in Fig. 10. 10b . Physically, the growth of dynamical heterogeneity with increasing 
q can be understood as the result of the competition between the length scale, q^ 1 , 
over which the dynamics is probed, and the typical particle displacement, <5, due to 
a rearrangement event. At very large q, qS > 1, one single event is sufficient to fully 
decorrelate g% — 1. In this regime, the instantaneous relaxation time depends on the 
time between successive events, which is a fluctuating quantity due to the Poissonian 
nature of the events. This yields very large fluctuations of the dynamics. As q de- 
creases, an increasing number of events is required to decorrelate gi — 1, thus leading 
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Fig. 0.10 Average dynamics and dynamical heterogeneity in a strong gel (adapted 
from (Duri and Cipcllctti, 2006]) with permission), a): q dependence of the relaxation time 
r r of the intensity correlation function (left axis, open circles) and of the stretching expo- 
nent (right axis, solid and semifilled squares). The data are normalized with respect to the 
parameters of the model (see also the text): S — 250 nm is the average particle displacement 
due to one single rearrangement event and 7 = (960 s) _1 is the event rate. Note that the 
motion is, on average, ballistic-like (r r ~ g -1 ) and that the decay of gi — 1 is steeper than 
exponential ("compressed" exponential, p > 1). The lines are the predictions of the model, b) 
Inset: dynamical susceptibility xili T ) f° r various q. Main plot: q dependence of the height, 
X* , of the peak of x( r )- The dashed line is a power law fit with an exponent 1.13 ± 0.11, the 
solid line is the model 

to smoother dynamics. 

Quite intriguingly, the main features of the average dynamics reported for the 
strongly attractive gels discussed above have been also found in a large variety of 



glassy soft materials (for a review, see e.g. Ref.(Cipelletti and Ramos, 20051). Mod- 
els inspired by Ref. ( |Duri and~C ipcllct ti, 2006] ) have been used also to describe the 
dynamics of nanoparticles embedded in molecular systems approaching the glass tran- 
sition ( |Caronna et ai, 2008| Guo et ai, 2009). Moreover, photon correlation imaging 



experiments (Sessoms et ai, 2009 Maccarrone et al., ) on several systems (from con 



ccntratcd soft spheres to repulsive glasses of charged clays and humidity-sensitive 
biofilms) suggest that system-spanning correlations of the dynamics may be a ubiqui- 
tous feature of jammed materials. This is in stark contrast with supercooled colloidal 
hard spheres, where spatial correlations of the dynamics extend over a few particle sizes 
at most, as discussed in Sec. 10.3.51 It is likely that the predominantly elastic behavior of 
jammed materials is responsible for the ultra-long spatial correlations of the dynamics 
observed in those systems, since the strain field set by a local rearrangement can prop- 
agate over large distances before being appreciably damped. Indeed, internal stress 
relaxation is often invoked as the origin of theses dynamics (Cipcllc tti et al., 2000| 
|Bouchaud and Pitard, 2002[ ), although a complete understanding of the physical mech- 
anisms underlying this peculiar yet general relaxation behavior is still lacking. 
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0.6 Perspectives and open problems 

While the average dynamics of glassy colloidal systems has been intensively studied 
since the 1980's, experiments on dynamical heterogeneities started only about twelve 
years ago, spurred by advances in microscopy and light scattering methods and stim- 
ulated by numerical works. A (partial) list of what we have learned in the past years 
on both the average dynamics and dynamical heterogeneity of glassy colloids includes: 

• Microscopy experiments have allowed us to observe directly what caged motion 
and cage rearrangements look like. Data have been used to quantify the meaning 
of "caging " QWeeks and Weitz, 2002] ). 

• Scattering experiments have shed new light on the dynamics of concentrated hard 
spheres, for which equilibrium dynamics above the ergodic-non ergodic transition 
predicted by mode coupling theory have been reported ( |Brambilla et ai, 2009[ ). 

• A variety of microscopy and light scattering experiments agree on both the ex- 
istence and the magnitude of dynamical heterogeneity in several colloidal systems, 



both 2D (Marcus et ai, 1999| Konig et al, 2005 Ebert et ai, 2009 Mazoyer et al, 2009) 



and 3D ( |Kasper et al., 1998[|Kegel and van Blaaderen, 2000||Weeks et al, 2000[), 

with both hard ( |Kegel and van Blaaderen, 2000| Weeks et al, 2000 ) and soft ( Konig et al, 2005[ 



|Yunker et ai, 2009 ; Scsso mFet al., 2009D repulsive interactions. Along with simu- 



lations dDoliwa and Heuer, 1998||Donati et al, 1998l|Glotzer, 2000[|Yamamoto and Onuki, 1998D , 



this provides nice evidence that dynamical heterogeneities are ubiquitous in glassy 
systems and not artifacts of one particular colloidal system, one particular exper- 
imental technique, or one particular simulation method. 

Experiments on colloidal hard spheres have provided some of the first experimen- 
tal quantitative evidence that the length scale of dynamical heterogeneities in- 
creases when approaching a glass transition (Weeks et at, 2000; Bcrthicr et ai. 2005; 



Weeks et al., 2007). 



• Light scattering experiments on deeply jammed attractive or soft repulsive sys- 
tems (Ballest a et ai, 20 08a; Ses soms et al., 20 09; Dur i et al., 20 09 b\ Ma ccarrone et al. 
have unveiled a richcr-than-cxpcctcd scenario, with ultra-long ranged spatial cor- 
relations of the dynamics not observed so far in simulations. 

In spite of these advances, several questions remain open, making dynamical het- 
erogeneity an exciting field of research: 

• Is there a structural origin of dynamical heterogeneity? While for diluted, at- 
tractive systems this has been shown to be the case, for concentrated, repulsive 
particles a clear answer is still lacking. Progress in this area will likely require 
the measurement of non-conventional structural quantities, e.g. the "point-to- 
set" correlation function introduced in numerical works ( |Biroli et al., 2008] ) and 
discussed in Chapter 2. 

• What is the behavior of dynamical heterogeneity in colloidal glasses and its rela- 
tionship with aging? While recent work on soft particles suggests that the slowing 
down of the dynamics during aging may be associated with a growth of spatial 



correlations of the dynamics (Yunker et al., 2009), this was not seen in experi- 



ments with harder spheres QCourtland and Weeks, 2003D ; this question requires 
further exploration. 
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• Is the non-monotonic tp dependence of dynamical fluctuations observed in some 
systems ( |Ballesta et al., 200 8b; Scs soms et al., 2009] ) a general feature? Although 

a somehow similar behavior has been reported for granular systems ( Lechenault et ai, 20081, 
the explanation proposed for colloids and grains arc different: can these contrast- 
ing views be reconciled? 

• Recent work ( [Duri et al, 2009 6[ |S essoms et al., 2009||Maccarrone et al., D on jammed 
soft materials suggests that the relatively short-ranged correlation of the dynam- 
ics of supercooled hard spheres may be the exception rather than the rule, since 
system-size dynamical correlations are observed in those materials. What is the 
physical origin of these correlations? Why do the numerical simulations not cap- 
ture these correlations? 

• Most work on colloidal systems has been devoted to hard sphere-like systems. A 
general understanding of the role of the interaction potential on the slow dynamics 
and dynamical heterogeneity is still lacking. Recent experiments show that softer 
colloids have strikingly different behaviors than hard colloids ( |Mattsson et ai, 2009 &D , 
although these results are not fully understood. 

Finally, we remark that numerical simulations can nowadays probe a range of relax- 
ation times comparable to that explored by experiments (see e.g. ( |Brambilla et ai, 2009D ). 
On the one hand, this calls for a more rigorous approach to the design of new exper- 
iments, since the question of what can be learned from experiments that simulations 
can not address needs to be asked. On the other hand, this opens the exciting possi- 
bility to compare in great detail numerical and experimental results, thereby allowing 
one to identify the physical mechanisms that are relevant in determining the slow 
relaxation of glassy colloidal systems. 
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